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0. Introduction 

A point process is an ensemble of random locally finite point configurations in 
a space X. A convenient tool for dealing with point processes is provided by the 
correlation functions (they are similar to the moments of random variables |19j). 
The correlation functions are indexed by the natural numbers n = 1,2,..., and the 
nth function depends on n variables ranging over X. A point process is said to be 
determinantal [5], [26] if its correlation functions of all orders can be represented 
as the principal minors of a correlation kernel K(x, y), a function on X x X. Then 
the whole information about the point process is contained in this single function 
in two variables. 

Concrete examples of determinantal point processes arise in random matrix the- 
ory, in some problems of representation theory, in combinatorial models of math- 
ematical physics, and in other domains. A remarkable fact is that the correlation 
kernels of determinantal point processes from various sources reveal a similar struc- 
ture. 

We will be interested in limit transitions for point processes depending on a 
parameter. For instance, a typical problem in random matrix theory consists in 
the study of limit properties of the spectra of random matrices of order N as the 
parameter TV goes to infinity. In the models of point processes arising in repre- 
sentation theory, the limit transitions may be related to the approximation of an 
infinite dimensional group by finite-dimensional ones [23] . [6]. 

The language of correlation functions is well adapted to studying limit transi- 
tions: the convergence of processes is controlled by the convergence of correlation 
functions (just as the convergence of random variables is controlled by the conver- 
gence of their moments). For determinantal processes, the situation is simplified, 
because the convergence of the correlation functions is ensured by the convergence 
of the correlation kernels. In the present note we will study just limit transitions 
for kernels. 

Certain correlation kernels K(x, y), which arise in many concrete examples from 
various domains, share the following common properties: 

• The space X, on which the kernel is defined, is a subset of R. 

• The kernel can be written in the so-called integrable form, which in the 
simplest variant looks as follows: 
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(0.1) K{x t y) = A M B M-_*to A <y), 

where A and B are certain functions on X. 

In such a situation, the convergence of kernels depending on a parameter is 
often extracted from the asymptotics of the functions A and B with respect to the 
parameter. This obvious way, however, may require hard computations if A and B 
are sufficiently complex special functions. □ In the present note, we discuss another 
approach, which exploits one more property of the kernels, which also holds in many 
concrete examples: 

• The operator K in the Hilbert space L 2 (X) (with respect to a natural 
measure), given by the kernel K(x,y), is a projection operator. Moreover, 
K can be realized as a spectral projection V(A) for a certain selfadjoint 
operator D acting in L 2 (X). Here A is a certain part of the spectrum of the 
operator D, and this operator is determined by a difference or differential 
operator on X. 

As soon as there is a natural link between K and D, the following simple idea 
arises: to deduce the convergence of the correlation kernels from that of the corre- 
sponding selfadjoint operators D. This idea was stated in 10J and then employed 
in [T3] and [3]. The aim of the present note is to demonstrate the efficiency of 
such an approach on a number of other examples. I consider only rather simple 
examples where the point processes live on the one-dimensional lattice and not on 
a continuous space, and the limit transition does not require a scaling. I think, 
however, that the method can be useful in more complex situations, too. Some 
indication on this is contained in §3.4 where it is explained how one can guess, in 
a very simple way, the scaling leading to the Airy kernel. 

The structure of the work is as follows. In §1 we recall general definitions related 
to point processes (for more detail, see pjj], [26 ). In §2 we introduce the main 
model: a 3-parameter family of probability measures on the Young diagrams 0). 
There it is also explained how to pass from measures on the Young diagrams to 
point processes on the lattice and where our main object, the "hypergeometric" 
difference operator on the lattice, comes from. In §3 we discuss limit transitions 
related to the degeneration of the z-measures to the Plancherel measure. In §4 
we study another limit regime that leads to the so-called Gamma kernel. This 
interesting kernel (it arose in [7]) describes the asymptotics of the fluctuations of 
the boundary of a random Young diagram near the point of intersection with the 
diagonal (in another language, limit properties of the lowest Frobenius coordinates). 

I am grateful to Alexei Borodin for discussions. I am also grateful to Rais 
Ismagilov and Yuri Ncrctin for valuable critical remarks. 



For instance, a solution of the problem of harmonic analysis on the infinite— dimensional 
unitary group [6] required computing the asymptotics of of Hahn type orthogonal polynomials, 
which are expressed through the hypergeometric series 3-F2- 

■^These measures, called the z— measures, were introduced in [5] for solving the problem of 
harmonic analysis on the infinite symmetric group. See also the survey paper [24] . 
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1. Definitions 

Let X be a topological space. By a configuration in X we mean a subset X a X 
without accumulation points. Q The space of all configurations in X is denoted as 
Conf(X). If X is discrete then Conf(X) is simply the set 2 X of all subsets of X. We 
will be dealing with probability Borel measures on Conf(X). Given such a measure 
P, one may speak about an ensemble of random configurations or a random point 
process. 

In what follows the space X will be discrete, so that all the next definitions will 
be given for this simple case. About the general case, see, e.g., [19]. 

We assign to P a sequence of functions on X, X x X, X x X x X, . . . , called 
the correlation functions. Here, by definition, the nth correlation function p^ is 
obtained in the following way: its value p^ n '(xi, . . . ,x n ) at an n-tuple of points 
equals if among these points there are repetitions; otherwise the value is equal 
to the probability of the event that the random configuration X contains all the 
points x%, . . . , x n . 

The initial measure P is determined by its correlation functions uniquely. 
A point process P is said to be determinantal [5] , [26] , if there exists a function 
K(x, y) on X x X such that for any n = 1, 2, . . . 

p (n) {xi, ...,x n ) = det[K(x ll x j )]l j=1 . 

Let us call K{x, y) the correlation kernel of the process P. The operator in 
the coordinate Hilbert space £ 2 (X) with the matrix [K(x,y)] will be called the 
correlation operator and denoted by K . 

The whole information on the determinantal point process P is contained in its 
correlation kernel (equivalently, correlation operator). 

A simple but important example of determinantal point processes is afforded by 
the orthogonal polynomial ensembles . Let us give the definition for the discrete case 
we are interested in (for more details, see the survey paper [18]). 

Let X C K be a discrete subset and W(x) > be a weight function defined on X. 
Fix N = 1,2,.... By definition, the ensemble in question is formed by the _/V-point 
configurations X = {xi, . . . , xn} C X, whose probabilities are given by the formula 

N 

(1.1) Prob(X) = const • JJ Wfa) ■ ][J (x l -x j ) 2 , 

i=l l<i<j<N 

where "const" is a suitable normalizing factor. 

It is well known (see, e. g., [18]) that such an ensemble is a determinantal process 
and its correlation kernel K(x,y) is expressed through the orthogonal polynomials 
Po = l,px,p2, . . . with the weight function W(x), as follows: 

K(x,y) = ^2 Pi{x)pi{y), x,y £ X, 

i=0 

where 

Pi(x) = y/W(x) , XGX, Z = 0,1,..., 

\m\\w 

4 More precisely, what we termed a configuration should be called a simple locally finite con- 
figuration. Here "simple" means that we exclude multiple points which are allowed in a more 
general definition. 
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and || • || w stands for the norm in the weighted Hilbert space £ 2 (X, W). Observe that 
the functions Pi{x) form an orthonormal system of vectors in £ 2 {X). Thus, the kernel 
K(x, y) is the matrix of a finite-dimensional projection in £ 2 (X): specifically, of the 
projection onto the subspace of polynomials of degree at most N — 1, multiplied by 

VW(xj. 

Since K(x, y) is equal to the product of y / W(x)W(y) with the iVth Christoffel- 
Darboux kernel for our system of orthogonal polynomials, the kernel K(x,y) can 
be written in the form 

= congt p N (x)p N ^(y)-p N ^(x)p N (y) 
x-y 

and hence can be represented in the integrable form (jO.ip . 

2. Z-MEASURES AND DISCRETE ORTHOGONAL POLYNOMIAL ENSEMBLES 

2.1. Z measures [5], [9]. Denote by Y the set of all partitions A = (Ai, A2, . . . ), 
which we identify with the corresponding Young diagrams. The z-measure with 
parameters 

zeC, z'eC, £eC\[i,+oo) 

is the (complex) measure M ZtZ >£ on the set Y assigning to a diagram A S Y the 
weight 

(2.1) M^(A) = (1 - 0"' e |A| (z)x(z'h (^r) 2 ■ 

Here |A| := ^ ■ A^ (cquivalcntly, |A| is the number of boxes in the diagram A); 

e(X) 

(2.2) (z) X = l[(z-i + l) Xi = H (z + j-i), 

i=i dJ)e\ 

where (a) m — a(a+ 1) ... (a+m— 1) is the Pochhammer symbol; £(X) is the number 
of nonzero coordinates A.;; the second product in (|2.2[) is taken over all boxes (i, j) 
of the diagram A, where i and j denote the numbers of the row and the column 
containing the box; finally, dim A denotes the number of standard tableaux of shape 
A (equivalently, the dimension of the irreducible representation of the symmetric 
group of degree |A| indexed by the diagram A). 
The following summation formula holds: 

Here we assume that either |£| < 1 (then the series absolutely converges) or the 
parameters z, z' are integers of opposite sign (then the series terminates). 

Note that the z-measures are a particular case of the Schur measures, see |22j). 

Note also two properties of the z-measures: 

(2.3) M 2 , Z ^(A) = M z , tZ , ( (\) M z . z ,. 6 (\ r ) = M_ z ,_y, ? (A), 

where A' denotes the transposed diagram. 

In what follows we consider only the z-measures with real nonnegative weights: 
Afz,z',£(A) > for all AeY. These are probability measures on Y. They fall into 
the following 4 series. 
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• Principal series: The parameters z and z' are complex conjugate and nonreal 
while the parameter £ is real and < £ < 1. 

• Complementary series: The both parameters z and z' are real and contained 
in one and the same open interval of the form (m, m + 1) where m £ Z, while the 
parameter £ is the same as above (0 < £ < 1). 

• First degenerate series: One of the parameters z, z' (say, z) is a nonzero integer, 
while the other parameter (then it will be z') is a real number of the same sign such 
that \z'\ > \z\ - 1. Again, < £ < 1. 

• Second degenerate series: The parameters z and 2/ are integers of opposite 
sign. The parameter £ is now subject to another condition: £ < 0. 

If a z-measure belongs to the principal or complementary series then (z + k) (z' + 
k) > for all integers k € Z, which implies that the weights of all diagrams are 
strictly positive. For the z-measures of the degenerate series some of the weights 
vanish but all nonzero weights are strictly positive. 

2.2. More details on degenerate z measures. The Meixner and Krawtchouk 
ensembles. Consider the first degenerate series. By virtue of the symmetry rela- 
tion (|2.3[) it suffices to take z = N and z' = N + c — 1, where N = 1,2,... and 
c > 0. Then the weight of a diagram A vanishes precisely when £(A) > N. Thus, 
the support of the z-measure is the set of the diagrams contained in the horizontal 
strip of width N. Such diagrams A are in a one-to-one correspondence with the 
A-point configurations L on Z + := {0, 1,2,...}: 

(2.4) A= (Ai,...,Aiv,0,0,...) < — ► L=(h,...,l N ), 

where 

h = Xi + N - i, 1 < i < A. 

It is readily checked that in terms of the correspondence (|2.4[) , the weight Mjy jv+ c _i ^(A) 
can be written in the form (jl.ip . where we have to set 2j = U and take as W the 
weight function for the Meixner polynomials |17j on the set X = Z + : 

j^Meixncr^ = ^ I £ Z+ . 

Thus, the correspondence A — > L converts the z-measure of the first degenerate 
series with parameters z = N, z' = N+c—1, £ into the JV-point Meixner orthogonal 
polynomial ensemble with parameters c and £. 

A similar fact holds for the second degenerate series. Let z = N and z' = —N' 
where N and N' are two positive integers, and assume £ < 0. Then the weight of a 
diagram A does not vanish if and only if £(\) < N and < N', that is, A has to 
be contained in the rectangular □jy.iV' with N rows and N' columns. The support 
of such a measure is finite. 

The same correspondence A «-> L as in (|2.4p gives a bijection between the dia- 
grams A C □jv,jv' an d the iV-point configurations L on the finite set {0, 1, ... , N}, 
where N = N + N' — 1. Then the z-measure of the second degenerate series turns 
into the A^-point orthogonal polynomial ensemble on X = {0,1,..., A} that is 
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determined by the weight function 



W KTawtchonk (l) = ()p l (l-p) N - 1 , l = 0,...,N, 
for the Krawtchouk orthogonal polynomials |17j with the parameters 
p:=^— , N:=N + N'-1, 





Note that the condition £ < guarantees that < p < 1 . 

Thus, the Meixner and Krawtchouk ensembles can be extracted from the z- 
measures as a rather particular case. As we will see, it is also helpful to use the 
inverse transition: 

Thesis 2.1. The z-measures can be obtained from the Meixner or Krawtchouk 
ensembles via analytic continuation (interpolation) with respect to the parameters. 

Let us explain this statement (for more detail, see [9] and [8]). Note that the 
support of a degenerate measure enlarges as the parameter N grows (or the two 
parameters N and N' grow) , so that any given diagram A £ Y can be covered. Fur- 
ther, as seen from (|2.1[) . for any fixed diagram A, its weight M z z ic(X) is an analytic 
function in £ whose Taylor coefficients at £ = are given by polynomial functions 
in z and z'. This allows one to extrapolate results about the degenerate series to 
the general case, because such functions in z, z' , and £ are uniquely determined by 
their restriction to the subset of the values of the parameters corresponding to the 
degenerate series. 

As we will see, Thesis 2.1 explains the origin of the difference operator (|2.6p 
introduced below. 

2.3. The hypergeometric deifference operator. Consider now the principal 
and complementary series. We cannot use the correspondence (|2.4|) anymore. In- 
stead of it we will introduce another correspondence, which will take arbitrary 
partitions A £ Y to semi-infinite point configurations X_ on the lattice U = Z + | 
of half-integers: 



(We call the configuration X_ semi-infinite because it contains all the lattice 
points that are sufficiently far to the left of zero and does not contain points suffi- 
ciently far to the right of zero.) 

There is a simple link between (12.4[) and (|2.5[) . Given a diagram A £ Y, let N be 
so large that N > £(X), so that the Appoint configuration L determined by (|2.4p 
exists. Then X_ is obtained from L by shifting it to the left by N — h and next 
adding the left-infinite "tail" {\ - i \ i = N + 1,N + 2, ...}. 

Let us assume the parameters (z, z') lie in the principal or complementary series. 
Consider the following difference operator on the lattice 11: 



(2.5) 



A 



X ={Ai-i+i | i= 1,2,...}. 



(2.6) V,s, t f{x) = ^{z + x + \){z' + x + \)f(x + 1) - {x + £(z + z' + x))f{x) 

+ ^(z + x-±)(z> + x-t)f(x - 1). 
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This difference operator was introduced in the papers [8] and [9] as the result of 
analytic continuation, in the spirit of Thesis |2~T1 of the Meixner difference operator. 
In more details, the analytic continuation procedure is as follows: 

1. Consider the Meixner difference operator on Z + : 

p Mci X ncr /(x) = fa + ^ + 1} _ (x + fa + c)) J {x) + x f {x _ 1} . 

It multiplies the nth Meixner polynomial by —(1 — see [17], formula (1.9.5). 
Add to pMoixncr the constant term (l — £)(iV — |) and observe that the first N 
Meixner polynomials are just those eigenvectors of the resulting operator 

(2.7) pMcixncr + (1 _ ^)(JV — i) 

that correspond to nonnegative eigenvalues. 

2. Pass from the weight £ 2 -space € 2 (Z + , W Mc,xncl ) to the ordinary space £ 2 (Z + ). 
This means that (|2.7|) should be replaced by the composition of operators 

(2.8) (^Meixner)! Q ^Meixner + ^ _ _ 1^ Q (^Meixner) - § _ 

3. Make in l|2.8[l the change of the variable I ^ x :— I — (N — |), which 
means that the support Z + of the weight function is transformed to the subset 
{ — (N — |) + m | m = 0, 1, 2, . . . } C U . Note that this subset grows together with 
N and in the limit it exhausts the whole lattice 11 . 

4. Formally substitute N — z and c — z' — z + 1 into the coefficients of the 
resulting operator. 

Then we get T> z ^ z i ^. Likewise, T> ZtZ '£ can be obtained from the difference oper- 
ator connected to the Krawtchouk polynomials. 

Denote by {e x }, x £ Z', the natural orthonormal basis in the Hubert space £ 2 (Z') 
and let ^q(Z') C ^ 2 (Z') stand for the algebraic subspace formed by linear combi- 
nations of the basis elements. We will consider T> Z;Z i^ (see ()2.6[) ) as a symmetric 
operator in i 2 (U) with domain £q(Z'). 

Proposition 2.2. The operator T> ' z iZ / ^ defined in this way is essentially selfadjoint. 

Proof. The standard way is to check that the eigenfunctions of the difference op- 
erator (|2.6p with nonreal eigenvalues do not belong to 1 2 {U). This can be done 
by computing explicitly the eigenfunctions (they are expressed through the Gauss 
hypergeometric function, see [9]). The computation can be simplified if one replaces 
(|2.6[) with the difference operator 

\/l\x+\\f(x + \) -(x + £(z + z' + x))f{x) + y/l\ x -\\f(x- 1), 

which differs from (|2.6[) by a bounded operator (adding a bounded selfadjoint op- 
erator preserves essential selfadjointness) . 

I will sketch now another argument. It is perhaps less elementary but it requires 
no computations and builds a bridge to representations of the group SL(2), which 
turns out to be useful in other questions related to the z-measures. 

Consider the complex Lie algebra gc = st(2, C) with the basis 



E = 



"0 


f 




"0 


0" 


, H = 


'1 " 





, F = 







1 





-1 
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and the commutation relations 

[H,E] = 2E, [H,P] = -2F, [E,F] = H. 

The real linear span of the elements E + F, i(E — F), and iH is a real form g C Qc- 
this is the Lie algebra su(l, 1). 

Define a representation S z ,z' of the Lie algebra gc in the algebraic space ^o(Z') C 
£ 2 (Z') by the following formulas for the basis elements 



S,,»'(E) e x = ^j( z + x+\){z' + x + \) e x+1 

(2 ' 9) S Z , Z ,(F) e x = - yf(z + x - \){z> + x - i) e x _! 

(H) = [z + z' + 2x)e x 

It is readily checked that the above operators satisfy the commutation relations and 
hence do determine a representation. 

The representation S z . z > equips with the structure of an irreducible Harish- 

Chandra module over q such that this algebra acts by skew symmetric opera- 
tors. It follows that S ZyZ ' gives rise to an irreducible unitary representation of 
the group ST/(1,1)~, the simply connected covering group of the matrix group 
SU(1, 1) = SL(2,M.) (see [T2], p. 14, where a proof is given for a concrete example 
but the argument holds in the general case). Moreover, according to a well-known 
theorem by Harish-Chandra, all vectors from £^(1') are analytic vectors of this 
representation. 

On the other hand, as seen from (12. 6|) and (12. 9| . 



(2.10) V z>z , t6 = S z , z , (v^E - V?F ^ i±i h) + i^(z + z')l, 

where 1 denotes the identity operator. Therefore, all vectors from ^q(Z') are ana- 
lytic vectors of the operator T> z . z i^, which implies ([25]. X.39) that it is essentially 
selfadjoint. □ 

Remark 2.3. The use of the representation S z g i of the Lie algebra s[(2,C) is 
suggested by Okounkov's work [21j . Note that the corresponding unitary represen- 
tation of the group SU(1, 1)~ belongs to the principal or complementary series, in 
exact agreement with the terminology for the z-measures (which was introduced 
prior to [21] . relying on a purely formal analogy). 

Let us agree to denote by -D z ,z',£ the selfadjoint operator in l 2 (U) that is the 
closure of the operator T> ZjZ i^. It follows from Proposition 12.21 that the domain of 
the operator D z , z >^ consists of all those functions / G £ 2 (Z') that remain in 1 2 {U) 
after application of the difference operator (|2.6|) . 



Proposition 2.4. The selfadjoint operator D z<z i^ has simple, pure discrete spec- 
trum filling the subset (1 — £)Z' C R. 

Proof. The matrix v^E — V? F — i£* H that enters (|2.10[) is conjugated with the 
matrix — ^^-H by means of an element of SU(X, 1). Therefore, D z z > ^ is unitarily 
equivalent to the closure of the operator S z ,z'{— ^^ll) + ^-^(z+z')l. This operator 
is the diagonal operator in the basis {e^} acting as multiplication of e x by — (1 — £)x. 
This makes the claim evident. □ 
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Thus, for any triple of parameters (z, z', £) from the principal or complementary 
series there exists an orthonormal basis {ip a -z,z',^} in ^ 2 (Z'), indexed by points 
a e U and such that 

^z,z'4^a;z,z'^ = (1 - QaipaWfr « € 

As shown in [9], the eigenvectors ipa;z,z',£ can be written down explicitly: for them 
there exists a convenient contour integral representation and also an expression 
through the Gauss hypergeometric function. 

In view of the connection of the functions ip a -z,z',^ with the hypergeometric func- 
tion I will call T> z z i ^ the hypergeometric difference operator. 

If A is a selfadjoint operator for which is not a point of discrete spectrum, then 
we will denote by Proj + (A) the spectral projection corresponding to the positive 
part of the spectrum of A. In particular, as the spectrum of D z ,z' £ does not contain 
0, we may form the projection Proj + (-D 2j z'.^). 

The image of the z-measure M ZyZ >£ under the correspondence A > X_ introduced 
in (|2.5p determines a point process on Z', which we denote as P z ,z< 

Theorem 2.5. P z ,z',£ * s a determinantal process. Its correlation kernel, denoted 
as K.z,z',i{x,y), is the matrix of the spectral projection Proj + (-D Z)2 ' ) {). That is, 

Kz,z',i( X >y) = ^ i>a;z,z',e{2)ll>a;z,z'Av)< X,yEZ'. 

(We do not put the bar over ip a ;z,z',z(y) because the eigenfunctions are real- 
valued, see formula (2.1) in [S].) 

In the above formulation, the result is contained in [9]. There are also references 
to previous works. We call K z>z / t ^(x,y) the discrete hypergeometric kernel. It can 
be written in the integrable form (|0.1[) . see [S], Prop. 3.10. 

3. PLANCHEREL MEASURE, DICRETE BESSEL KERNEL, AND DISCRETE SINE 

KERNEL 

3.1. Poissonized Plancherel measure. Consider the limit regime 
(3.1) £ — » 0, z — > ao, z — > oo, £zz' — ► 8, 

where 8 > is a new parameter. Then the z-measures M ZiZ i t £ converge to a 
probability measure Mg on Y: 

M,(A)= ii imM„„, £ (A) = e-»^ (^) 2 . 

The measure Mg arises in the result of poissonization of the sequence the Plancherel 
measures JlfW, n= 1,2,..., where 

M{n) (dimA) a |A| = n _ 
n! 

For this reason Mg is called the poissonized Plancherel measure with parameter 0. 

For more details about the Plancherel measures JlfW see the papers |20j and 
[27] : see also [16] and [H]. The poissonized version was first considered in [I]; 
further results were obtained in [4] and [IS] . 
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As in §2, we use the correspondence A i— ► X_ to pass from the measure Mg on Y 
to a point process on 11 \ let us denote this point process as Pg. It is determinantal; 
its correlation kernel is described in §3.2 below. 

3.2. From discrete hypergeometric kernel to discrete Bessel kernel. In 
the regime (|3.1[) . the formal limit of the difference operator T> z ,z',$ is the difference 
operator X>f csscl on the lattice U , acting according to formula 

vf essel f(x) = Vef(x + 1) - xf(x) + Vef(x - 1). 

Proposition 3.1. Consider Dg essel as a symmetric operator in £ 2 (Z') with domain 
£q(Z'). Then this is an essentially selfadjoint operator. 

Proof. Indeed, Vf CSSQl is the sum of the diagonal operator / t— ► —xf and a bounded 
operator. □ 

Alike the hypergeometric operator (|2.6|) . the operator X>f cssel is connected with a 
group representation, only the group is different, namely, it is the universal covering 
group G for the group G of motions of the plane K 2 . Let us discuss this connection 
shortly. 

Let g be the Lie algebra of the group G and gc be the complexification of g. In 
0C, there is a basis {E, F, H} with the commutation relations 

[H,E]=E, [H,F] = -F, [E,F]=0. 

Consider the representation S of the Lie algebra gc in the dense subspace ^o(Z') C 
£ 2 (Z') defined on the basis elements by the formulas 

S(E)e x = e x+ i, 5(F)e x = -e x _i, S(H)e x = xe x , 

It is directly checked that all vectors from ^g(Z') are analytic vectors for S and g acts 
by skew symmetric operators. Therefore, S gives rise to a unitary representation 
of the group G. We have 

which is an analog of equality (|2 . 1 0|) . 

Denote by Df essel the closure of the symmetric operator T>f csscl . According to 
Proposition 13. 11 Dg is selfadjoint. 

Set 

i>a-,e{x) := J x+a (2Ve), xez', 

where a 6 U is a parameter and J m ( • ) is the Bessel function of index m. 

Proposition 3.2. For any fixed 6 > 0, the functions ip a; g form an orthonormal 
basis in £ 2 (Z') and 

D 3e3sel i> a . t g = atp a .,g, a£Z'. 
Thus, jjBessei ^ as p Ure discrete, simple spectrum filling the lattice Z'ctt. 

Proof. The element v^E— v^F — H S ig is conjugated with — H under the adjoint 
action of the group G. It follows that the selfadjoint operators S(V0~E — \fd F — H) 
and S{— H) (where the bar means closure) are conjugated by means of a unitary 
operator corresponding to an element of the group G. Since D^essei coincides with 
S^v^E — vpF — H), its spectrum is the same as that of the diagonal operator 
S(— H). Thus, we see that the spectrum is indeed the same as was stated. 
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The fact that the functions ip a ;e(x) are eigenfunctions of the operator X>f csscl 
follows from the well-known recurrence relations for the Bessel functions ([2], 7. 2. 8 
(56)). Finally, the relation 1 1 V^o. ; 6> 1 1 2 = 1 can be proved in the same way as Propo- 
sition 2.4 in [9], by using the contour integral representation of the Bessel func- 
tions. □ 



Proposition 3.3. In the regime (|3.ip . the operators D z , z ' converge to the operator 
Df csscl in the strong resolvent sense. 

About the notion of strong resolvent convergence see the textbook [25], Section 
VIII.7. 

Proof. Obviously, D z ^ — > £)^ csscl on the subspace £q(Z'). By virtue of Proposi- 
tions [2T2] and [23 ^o(^') i s a common essential domain for all selfadjoint operators 
under consideration. Then the claim follows from a well-known general theorem 
([25], Theorem VIII.25). □ 

Corollary 3.4. The spectral projection operators Proj + (Z? 2iZ '^) strongly converge 
to the spectral projection operator Proj + (Z?g csscl ). 



Proof. This is a direct consequence of Proposition 13 . 31 and another general theorem 
([25], Theorem VIII. 24 (b)). For applicability of that theorem it is important that 
not be a point of discrete spectrum of the operators D zz i^ and Df osscl , and this 
fact follows from the description of their spectra, see Propositions 12.41 and 13.21 □ 

Corollary [3U shows that Proj + (L>f csscl ) serves as the correlation operator for 
Pq. The kernel of the projection Proj + (Z3f csscl ) can be written in the form 

*P-W) = £ = vg*-*^ w -*<* ( "*-^» 

— ' x — y 

aeZ' + 

This kernel is called the discrete Bessel kernel; it was independently derived in [15] 
and (in somewhat different form) in [?]. 

3.3. From discrete Bessel kernel to discrete sine kernel. Consider the fol- 
lowing limit regime depending on the parameter c G (—1, 1): 



(3.2) 6^oo, x^2cVe + x, xeZ. 

This means that as 9 goes to oo, we shift the lattice 11 so that to focus on a 
neighborhood of the point 2cV0. The following formal limit holds 

with the following difference operator in the right-hand side: 

p sinc /( ~ ) = /( ~ + 1} _ 2c/( ^ + /( ~ _ 1); ~ G z _ 

Proposition 3.5. Consider T> s c lnc as a symmetric operator in £ 2 (Z) with domain 
£q(Z). Then this is an essentially selfadjoint operator. 

Its closure, denoted as D^ lnc , has pure discrete spectrum of multiplicity 2, filling 
the interval (—2 — 2c, 2 — 2c). 
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Proof. The difference operator D^ ne is invariant under shifts of 1, which makes its 
study quite simple. Indeed, pass from the lattice Z to the unit circle |£| = 1 by 
means of the Fourier transform. Then our operator will become the operator of 
multiplication by the function 

C - 2(RC-c), 

which easily implies the claim about the spectrum. □ 

Proposition 3.6. In the regime (|3.2p , the operators £)f csscl converge to the oper- 
ator D s c mQ in the strong resolvent sense. 

The proof is the same as in Proposition 13.31 As in Corollary 13.41 we deduce 
from this the strong convergence Proj,(Dg csscl ) — > Proj + (£>^ mc ). The kernel of 
the projection Proj,(.D;! me ) is readily computed: it has the form 

= sin((arccosc)(,-^ 
n{x-y) 

and is called the discrete sine kernel [4]. 

3.4. From discrete Bessel kernel to Airy kernel. The Airy kernel first emerged 
in random matrix theory: under suitable conditions on the matrix ensemble, this 
kernel describes the asymptotics of the spectrum "at the edge", as the order of 
matrices goes to infinity. The Airy kernel is defined on the whole real line and is 
expressed through the classical Airy function Ai(u) and its derivative: 

rWUrw \ f + °° ... v.., ^ Ai(u)Ai'(v)-Ai'( U )Ai(v) 

K A " y (u,v)= / Ai(u+s)Ai(v+s)ds = — — — — — — , 

Jo u - v 

As shown in [4] and [15] , the discrete Bessel kernel Kf essel (x, y) converges to the 
Airy kernel K Any (u : v) if 9 — > +oo and the variable x € 11 is related to the variable 
u £ K by the scaling x = 26 1 / 2 + • u. Below we present a heuristic derivation 
of this claim. 

Let g(u) be a smooth function on M. Assign to it the function f{x) on 11 by 
setting fix) = g(u), with the understanding that x and u are related to each other 
as above. Then we get 

f{x ± 1) = g(u ± e^/ Q ) « giu) ± e^ 6 ■ g'iu) + \ 9^ ■ g"(u), 

so that 

Vf csscl f(x) = 9 1 ' 2 {f{x + 1) - x9-^f{x) + f(x - 1)) « 9^ig"iu) - ug(u)). 

This simple computation explains the origin of the factor 2 and the exponents 1/2 
and 1/6. 

Next, the eigenvalue equation 

turns, after the renormalization a = 9 1 ^ e s : into the equation 

V Miy iP = atp, se R, 

where 

V A ^giu) = g"iu)- U giu) 
is the Airy operator. Its eigenfunctions have the form 

%l> a {v) = Ai(u + s), set, 
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and the Airy kernel corresponds to the spectral projection onto the positive part 
of the spectrum. 

It would be interesting to make this argument rigorous by using the notion of 
Mosco convergence. 

4. From discrete hypergeometric kernel to Gamma kernel 

Let us fix the parameters {z, z') from the principal or complementary series and 
let £ go to 1. As shown in [7], in this regime, there exists a limit for the correlation 
kernels. Below we prove this fact in a new way (Corollary 14.41) . 

The convergence of the kernels implies that the images of the measures M z>z r £ 
under the correspondence (|2.5[) converge to a probability measure on Conf(Z') 
generating a determinantal point process. Note that for the initial z-measures on 
Y the picture is different: they converge to 0. Of course, this cannot happen on the 
compact space Conf(Z') = 2 Z . 

The first indication on the existence of a limit is the fact that the hypergeometric 
difference operator (|2.6p still makes a sense under the formal substitution £ = 1: 
then we get the difference operator 



(4.1) VI— f{x) = ^{z + x + \)(z' + x + \)f{x + l)-(z + z> + 2x)f(x) 

+ ^(z + .T-i)(z'+X-i)/(.T - 1). 

In the sequel we follow the scheme we have already worked out. 

Proposition 4.1. Consider D^ a ™ a as a symmetric operator in ^ 2 (Z') with dome 
£q(Z,'). Then this is an essentially selfadjoint operator. 



Proof. The both arguments given in the proof of Proposition 12.21 hold without 
changes. □ 

Denote by jj)S a ™ ma the selfadjoint operator that is the closure of the operator 

/rjgamma 

Proposition 4.2. £)8 a ™ ma /j as p Ure continuous spectrum. 

Proof. Note that for £ = 1, the relation (|2.10p turns into 

V &— = Sz z , ( E - F - H) . 

Observe now that E — F — H is a nilpotent element in isu(l, 1) C st(2, C). For any 
r > the matrix r(E — F — H) is conjugated to E — F — H by means of an clement 
from SU(1, 1). Therefore, r £)| a " lma is unitarily equivalent to rj)S a ™ ms \ This implies 
the claim. □ 

Actually, one can say more: the spectrum is simple and it fills the whole real 
axis. This can be deduced, for instance, from a suitable realization of the represen- 
tation S ZtZ i. The eigenfunctions of the operator £>s a ™ ma constituting the continuous 
spectrum can be written down explicitly, see [8]; they are expressed through the 
classical Whittaker function. 

Proposition 4.3. As £ — » 1, the selfadjoint operators D ZjZ /£ converge to the self- 
adjoint operator £>s a ™ ma j n ^ e strong resolvent sense. 
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The proof is the same as in Proposition 13.31 Just as in Corollary [33] we deduce 
from this 

Corollary 4.4. As £ — > 1, the projection operators Proj _l_(D ZjZ /^) strongly converge 
to the projection operator Proj + (Z?^ a ™ ma ). 

The kernel of the limit projection operator is explicitly described in [7j, [5]. 
It can be written in the integrable form (jO.ip . where the functions A and B are 
expressed through the Gamma function, which explains the origin of the term 
"Gamma kernel" . 
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